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ABSTRACT 

We report the discovery of PSR J2022+3842, a 24 ms radio and X-ray pulsar in the supernova 
remnant G76.9+1.0, in observations with the Chandra X-ray telescope, the Robert C. Byrd Green 
Bank Radio Telescope, and the Rossi X-ray Timing Explorer ( RXTE ). The pulsar’s spin-down rate 
implies a rotation-powered luminosity E = 1.2 x 10 38 erg s^ 1 , a surface dipole magnetic field strength 
B s = 1.0 x 10 12 G, and a characteristic age of 8.9 kyr. PSR J2022+3842 is thus the second-most ener- 
getic Galactic pulsar known, after the Crab, as well as the most rapidly-rotating young, radio-bright 
pulsar known. The radio pulsations are highly dispersed and broadened by interstellar scattering, and 
we find that a large ( Sf/f « 1.9 x 10~ 6 ) spin glitch must have occurred between our discovery and 
confirmation observations. The X-ray pulses are narrow (0.06 cycles FWHM) and visible up to 20 
keV, consistent with magnetospheric emission from a rotation-powered pulsar. The Chandra X-ray 
image identifies the pulsar with a hard, unresolved source at the midpoint of the double-lobed radio 
morphology of G76.9+1.0 and embedded within faint, compact X-ray nebulosity. The spatial relation- 
ship of the X-ray and radio emissions is remarkably similar to extended structure seen around the Vela 
pulsar. The combined Chandra and RXTE pulsar spectrum is well-fitted by an absorbed power-law 
model with column density Nr = (1.7 db 0.3) x 10 22 cm -2 and photon index T = 1.0 ± 0.2; it implies 
that the Chandra point-source flux is virtually 100% pulsed. For a distance of 10 kpc, the X-ray 
luminosity of PSR J2022+3842 is Lx(2-10 keV) = 7.0 x 10 33 erg s -1 . Despite being extraordinarily 
energetic, PSR J2022+3842 lacks a bright X-ray wind nebula and has an unusually low conversion 
efficiency of spin-down power to X-ray luminosity, Lx/E = 5.9 x 10 -5 . 

Subject headings: pulsars: individual (CXOU J202221.68+384214.8, PSR J2022+3842) — stars: neu- 
tron — supernova remnants — X-rays: stars 


1. INTRODUCTION 

The past decade has been witness to a highly produc- 
tive back-and-forth relationship between X-ray and radio 
studies of supernova remnants (SNRs). Their nonther- 
mal radio emissions, serving as tracers of particle acceler- 
ation, have motivated follow-up X-ray observations, most 
notably with the Chandra telescope. In many instances, 
X-ray imaging has revealed point-like sources surrounded 
by extended synchrotron structures, typically axisym- 
metric; these pulsar wind nebulae (PWNe) arise from 
the shocked outflows of relativistic particle winds driven 
by a central engine, the powerful time-varying fields in a 
neutron star magnetosphere. Targeted radio periodicity 
searches of PWNe have then uncovered previously un- 
seen pulsations, providing the spin periods and, eventu- 
ally, estimated ages, magnetic field strengths, and spin- 
down luminosities of the neutron stars, without which 
a meaningful physical understanding of the surrounding 
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phenomena would be beyond reach. A few specific ex- 
amples of this general picture of radio/X-ray symbiosis 
have been the discoveries of the pulsars in the remnants 
3C58, G54.1+0.3, G292.0+1.8, and in the “Mouse” (see, 
e.g., Camilo 2004). 

The Crab Nebula has been viewed for decades as a 
prototype, but there is growing evidence of PWNe with 
substantially different properties. At least three objects, 
among them the Vela SNR, have nonthermal radio spec- 
tra that are steep in comparison to the Crab (i.e., with 
spectral indices a & 0.6, where flux S v oc v~ a at fre- 
quency u) and morphologies characterized by radio emis- 
sion lobes straddling a central depression. While the na- 
ture of one such candidate PWN, DA 495, was once in 
question, X-ray imaging observations led to the discovery 
of an X-ray nebula and a presumed neutron star, both 
residing near the radio “hole” (Arzoumanian et al. 2008; 
hereafter, ASL+08). Radio and X-ray pulsation searches 
have not, however, detected a pulsar signature, so a key 
ingredient is lacking from our understanding of DA 495. 
In this paper, we describe a similar series of observations 
directed toward understanding the third member of this 
group of unusual SNRs, G76.9+1.0. 

SNR G76. 9+1.0 was discovered in a 408 MHz DRAO 
survev of the Cygnus X region and initially characterized 
as a “possible galactic object” (Wendker et al. 1991). 
Follow-up three-band radio observations using the Very 
Large Array (VLA) by Landecker et al. (1993) resolved 
a two-iobed structure 4' in size, with a connecting bridge 
of emission, embedded in a faint, roughly circular emit- 
ting region 9' x 12' across. Based on its steep, nonthermal 
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Figure 1. Left: Chandra ACIS-S3 0.5-7 keV X-ray image of the center of G 76.9+1.0, smoothed, exposure-corrected, and scaled to 
emphasize the field point sources. The 1/ x 1/ box is centered on CXOU J 20222 1.68 +3 842 14. 8, while the circled sources were used to 
register the image. Lower inset: The unsmoothed count-rate map within the box containing the candidate neutron star and PWN. 
Upper inset: The result of adaptive smoothing of the same image. Overplotted are the regions used to extract spectral information 
for the unresolved source (circle) and diffuse emission (ellipse). Right: 1.49 GHz map of SNR G76.9+1.0 obtained using the VLA (see 
Landecker et al. 1993 for details) with the X-ray pulsar and nebula regions superposed (contours). The pulsar lies along the ridge of 
emission connecting the two lobes, adjacent to the central “hole.” The arrangement is strikingly similar to the radio/X-ray configuration 
of the Vela pulsar and PWN (see Figure 10 of Dodson et al. 2003) as well as of DA 495 (Figure 1 of ASL-j-08). 


spectrum, polarization, and morphological similarities to 
center-filled SNRs (in important respects such as the lack 
of an c uter boundary and radially decreasing flux away 
from the central depression), these authors favored an 
SNR interpretation. In this picture, the radio lobes trace 
the Crab-like synchrotron emission of an evolved pulsar 
wind nebula. 

Here we present imaging, spectral, and timing studies 
of G76.9+1.0 in X-rays and radio, culminating in the dis- 
covery of PSR J2022+3842 and of its associated X-ray 
PWN. In §2, we describe analysis of a Chandra obser- 
vation that yielded the first identification of the SNR’s 
central engine, a neutron star with its surrounding non- 
thermd nebula. In §3, we describe the discovery of the 
radio pulsar J 2022+3842 and follow-up radio timing ob- 
servations with the Green Bank Telescope (GBT) that 
have provided its spin-down rate and pulse properties. In 
§4, we present the results of a deep Rossi X-ray Timing 
Explorer ( RXTE ) observation: we show that the pulsed 
X-ray flux is consistent with the Chandra point-source 
flux, cementing the pulsar’s association with the SNR. 

For reasons developed below, we adopt a distance to 
G76.9+ 1.0 of 10 kpc throughout this work. 

2. Chandra OBSERVATIONS AND RESULTS 

SNR G76.9+1.0 was observed for 54 ks by the Chandra 
observatory on 2005 August 01 UT using the Advanced 
CCD Imaging Spectrometer (ACIS; Burke et al. 1997) 
operating in the full-frame TIMED/VFAINT exposure 
mode (ObsID #5586). This detector is sensitive to X- 
rays in the 0.3-12 keV energy range with a resolution of 
A E/E ~ 0.06 FWHM at 1 keV. The CCD pixel scale is 


0.5", comparable to the telescope’s on-axis spatial reso- 
lution. Low-level data were reprocessed with the latest 
calibrations using the CIAO script chandra_repro. The 
data analysis made use of the CIAO (V4.3), FTOOLS 
(V6.10), CALDB (V4.4.2), and XSPEC (V12.6.0q) soft- 
ware packages. A total of 53.7 ks of live- time was ac- 
cumulated. No additional time filtering was necessary 
as the background rate was stable over the course of 
the observation. We followed the CIAO online science 
threads to create ACIS spectra and images. To al- 
low for imaging with the best available angular resolu- 
tion, we used the EDSER reprocessing option, applying 
an energy-dependent sub-pixel event-repositioning algo- 
rithm to each photon (Li et al. 2004). 

The left-hand panel of Figure 1 presents the exposure- 
corrected image, in the 0.5-7 keV band, acquired by the 
ACIS-S3 and S2 CCDs, which together covered the en- 
tire extent of the radio SNR. No morphological evidence 
of G76.9-K.0 itself is apparent nor is there any hint of 
resolved diffuse emission in images created in the soft 
(< 2 keV) or hard (> 2 keV) energy bands; instead, 
several faint point-like sources are detected. The bright- 
est of these, enclosed within a box in the figure, is lo- 
cated close to the optical axis and contains N = 1190 
photons. This source is not associated with any known 
stellar counterpart or object in the 2MASS image of the 
field. Its identification as a possible neutron star associ- 
ated with G76. 94-1.0 is reinforced by the presence of faint 
but distinct X-ray nebulosity surrounding a brighter, un- 
resolved core. To highlight the diffuse emission, we show 
as insets in Figure 1 both the raw and smoothed 0.5-7 
keV count-rate maps for the boxed 1' x 1' region cen- 
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tered on the source. The smoothed image was produced 
using a circular convolution kernel of varying size, requir- 
ing a minimum of 5 counts within the kernel diameter. 
Although the putative nebula is faint, there is a clear 
excess of emission elongated on a line with position an- 
gle ~ 26° North through East. This and the relative 
faintness of the point source argue against identification 
of the diffuse emission as a dust-scattering halo; these 
are normally circular and found around much brighter 
sources. 

The 1.49 GHz Very Large Array (VLA) radio image 
shown in the right-hand panel of Figure 1 places the 
X-ray nebula in the context of the double-lobed radio 
structure. The unambiguous alignment of the long axis 
of the diffuse X-ray emission with the radio ridge con- 
necting the lobes suggests that the radiation processes 
underlying the emissions in both bands are related, such 
as in the canonical picture of synchrotron emission from 
relativistic particles originating in the pulsar’s magneto- 
sphere. 

To refine the location of the putative pulsar, the 
X-ray image was registered using the coordinates of 
eight 2MASS infrared counterparts obtained from the 
NASA/IPAC Infrared Science Archive. The updated 
centroid is located at (J2000.0) R.A. = 20 h 22 m 2T689, 
Deck = +38 0 42'14782 with a la uncertainty of 
(0709,0707) in the two coordinates, respectively. The 
required shift in R.A. and Deck was (—07012, -07162). 
This source, CXOU J20222 1.68+3842 14. 8, is hereafter 
referred to by the shortened name CXOU J202221. 

We extracted a spectrum for CXOU J202221, and 
generated an appropriate response matrix, with the 
specextract CIAO script, using a circular 275-radius 
aperture (representing > 90% encircled energy, depend- 
ing on photon energy) centered on the source. The dif- 
fuse emission and cosmic and detector backgrounds con- 
tribute negligibly (~ 10 counts) in this region. With 
a maximum count rate in a pixel of < 4 x 10 -3 s -1 , 
photon pile-up could safely be ignored. The 0.5-10 keV 
spectrum was grouped with a minimum of 15 counts 
per spectral channel and fitted using XSPEC to an ab- 
sorbed power- law model, a common property of young, 
rotation-powered pulsars. (Blackbody models are for- 
mally allowed, but produce unnaturally high temper- 
atures.) An excellent fit, with x 2 — 09.5 for 70 de- 
grees of freedom (DoF) was obtained for an absorbing 
column Nn = (1.6 =b 0.3) x 10 22 cm“ 2 and photon in- 
dex T = 1.0 ± 0.2 (90% confidence intervals are used 
throughout). The absorbed 2-10 keV source flux is 
FpgR = 5.3 x 10“ 13 erg cm -2 s -1 , implying an isotropic 
luminosity of Lx = 7.0 x 10 33 erg s -1 at 10 kpc distance. 
As described in §4.2, this Chandra spectrum was also 
used in joint fits with RXTE data. 

To investigate the spectrum of the nebula, we used an 
elliptical extraction region (Figure 1), with semi-major 
and -minor axes of length 873 x 571, centered on the 
pulsar, wdth the point-source region ( r < 275) excluded. 
For the local background, we extracted photons from a 
concentric annular region well outside the nebular extent 
(20" < r < 40"). We also performed a MARX simulation to 
determine the small (2.6%) flux contribution within the 
nebular region from point-source photons scattered into 
the point-spread function wings; this simulation took 


into account the measured pulsar spectrum. The tally 
of background-subtracted counts from the nebula is then 
N = 76 ± 13. The faintness of the nebula does not allow a 
well-constrained spectral fit; instead we fixed the column 
density to the best-fit value for the unresolved source 
and assumed a nominal absorbed power-law model with 
T = 1.4, a value derived from the empirical trend relating 
the spectral indices of energetic pulsars and their nebu- 
lae (Gotthelf 2004). The absorbed 2-10 keV flux for the 
putative PWN is Fpwn » 4 x 10 -14 erg cm~ 2 s” 1 . The 
flux ratio Fpwn/Fpsr ~ 0-08 is well below that expected 
for similarly energetic pulsars (Gotthelf 2004). The un- 
certainty in this result is dominated by the small number 
of nebula photons and not the choice of T. 

3. RADIO OBSERVATIONS AND RESULTS 

Motivated by the Chandra discovery of the point 
source and possible PWN at the heart of SNR 
G76.9+1.0, we searched for radio pulsations with the 
GBT. The Spectrometer SPIGOT data-acquisition back- 
end (Kaplan et ah 2005) was used for 5 hours on each of 
two days, with 600 MHz of useful bandwidth centered 
at 1.95 GHz. The instrument provided 768 channels 
with 16-bit sampling of summed polarizations at 81.92 /is 
intervals. The PRESTO software (Ransom et al. 2002) 
was used to search for pulsations, and a candidate 24-ms 
periodicity was detected on both days with formal sig- 
nificances of 7a on MJD 54397 and 12 a on MJD 54400. 
Confirmation observations were carried out on 2009 May 
6 (MJD 54957) using the GUPPI data-acquisition sys- 
tem configured similarly to the discovery instrumentar 
tion, but with the addition of full-Stokes sensitivity. The 
difference in the pulse periods derived at these epochs 
gave the first indication of the pulsar’s spin-down rate 
and thence its canonical age, magnetic field strength, and 
spin-down power (see below). A search for single “giant” 
pulses yielded no compelling candidates, despite the fact 
that the pulsar’s estimated magnetic field strength at the 
light cylinder, suggested to be relevant in the production 
of giant pulses (Cognard et ah 1996), is approximately 
70% that of both the Crab pulsar and PSR B1937-I-21, 
and comparable to that of PSR B1821-24, all of which 
exhibit giant pulses. 

Additional radio observations, also with the GUPPI 
system, were carried out at roughly 3-month intervals 
to constrain the pulsar’s long-term timing behavior. 
Notably, the pulse periods measured during and after 
the confirmation observation differed substantially from 
those detected during the two discovery observations 
made 1.5 years earlier, implying that a spin “glitch” of 
magnitude A P/P ~ 1.9 x 10“ 6 had occurred at an un- 
known epoch in this interval. The spacing of the ra- 
dio timing observations together with apparent “tim- 
ing noise” instability in the pulsar’s rotation on simi- 
lar timescales precluded derivation of a phase-connected 
timing solution. We therefore determined the long-term 
average spin-down rate, 4.3192(27) x 10“ 14 , through a 
least-squares fit to the multi-epoch period measurements 
(Figure 2), a result consistent with the short-term X-ray- 
derived ephemeris shown in Table 1 (see §4.1). 

The pulsar’s flux density at 2 GHz, S 2 GHz> is 
60 its 4.8 GHz flux density of « 45 /iJy sug- 
gests that PSR J2022+3842 has an unusually flat spec- 
trum, a ~ -0.33. The implied radio pseudo-luminosity 
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Figure 2. Top: Radio pulse properties of PSR J2022+3842 aver- 
aged over a 700 MHz band centered at 1.95 GHz, and representing 
a total integration time of approximately 24 hours. The evolu- 
tion with pulse phase of total ( black line), linearly polarized {red), 
and circularly polarized {blue) flux are shown. Phase zero is ar- 
bitrary. Bottom: Radio pulse period evolution {red squares) with 
lcr uncertainty error bars; the blue cross corresponds to the RXTE 
observation (uncertainty much smaller than the plotted symbol). 
Residuals are relative to the timing model shown in Table 1 . The 
pulse period at the discovery epoch of the radio pulsations, around 
MJD 54400, was significantly longer than at later times, indicating 
that a spin-up glitch occurred at some point before the confirma- 
tion observation near MJD 54950. 

Si^GHzd 2 — 7D 2 0 HiJy kpc 2 is unremarkable, low in 
comparison to the majority of known pulsars, but an or- 
der of magnitude higher than the luminosities of other 
young, faint pulsars discovered in deep searches of SNRs 
(Camilo 2004). Even with a long cumulative exposure, 
our best average pulse profile (Figure 2) retains signifi- 
cant statistical noise; despite this faintness, a number of 
relevant pulse shape and polarization results are avail- 
able. 

The radio pulsations of PSR J2022+3842 are signifi- 
cantly affected by dispersion and scattering due to prop- 
agation in the interstellar medium, which are typically 
important at frequencies < 1 GHz, even in our 2 GHz 
observations. To estimate dispersion measure (DM) and 
the pulse broadening timescale r scatt , we formed a pulse- 
shape model assuming a gaussian-shaped intrinsic pro- 
file, a one-sided exponential-decay impulse response func- 
tion for the scatter- broadening, and a frequency de- 
pendence for the width of the exponential. We fit the 


S-band (i.e., 2 GHz) data to this model in four evenly- 
spaced frequency sub-bands, allowing the width of the 
intrinsic gaussian profile to vary and also accounting for a 
DM bias caused by the frequency-dependent pulse broad- 
ening due to scattering. This simple model fit the data 
well, but we caution that the assumptions inherent in the 
model may introduce significant systematic errors, per- 
haps as large as twice the statistical uncertainties quoted 
below. 

The NE2001 Galactic electron distribution model 
(Cordes & Lazio 2002) fails to accommodate the de- 
rived DM = 429.1 ± 0.5 pc cm -3 in the direction of 
PSR J2022+3842— formally, the model suggests a lower 
bound on the distance to the pulsar of 50 kpc. Instead, 
a likely over-density of free electrons in the Cygnus re- 
gion, along the line of sight, accounts for the higher-than- 
expected dispersion — a similar explanation has been ad- 
vanced for the unexpectedly large DM of the nearby PSR 
J20214-3651 (Roberts et al. 2002). 

In this direction, the Perseus Arm lies at a distance of 
approximately 6 kpc, and the Outer Arm at > 10 kpc. 
For convenience in scaling distance-dependent parame- 
ters and because the larger distance brings X-ray efficien- 
cies more approximately in line with those of other young 
pulsars (see §5), we adopt a distance d = 10 Dio kpc. In- 
dependent supporting evidence for the large distance de- 
rives from H I absorption and X-ray studies of the P WN 
CTB 87 (Kothes et al. 2003), roughly 2° from G76. 9+1.0 
on the sky. At a distance of 6.1 kpc, CTB 87 lies in the 
Perseus arm and its X-ray-derived absorbing column is 
Nil = (1.4 ± 0.2) x 10 22 cm -2 at 90% confidence (Safi- 
Harb, Matheson & Kothes, in preparation); although 
their uncertainties are large, the nominal Hi column for 
G76.9+1.0 is somewhat larger than that for CTB 87. A 
radio bright extragalactic point source adjacent to CTB 
87 exhibits a deep absorption component at a velocity of 
about —85 km s -1 not seen in the absorption profile of 
CTB 87. This feature could account for the larger fore- 
ground absorption toward G76.9+1.0 if the latter is lo- 
cated at a significantly larger distance, such as the Outer 
Arm or beyond it. 

The NE2001 model also predicts a timescale for pulse 
broadening due to interstellar scattering of 0.3 ms at 1 
GHz; instead, we find = 55 ± 7 ms, which scales 

to 3.8 ms at the center of our observing band but varies 
strongly across it. The long decay of the pulse, here 
averaged over the 700 MHz band, can be seen in the 
upper half of Figure 2. Consistent with this result, an 
exploratory observation centered at 800 MHz sky fre- 
quency failed to detect the pulsar; similar attempts at 1.4 
GHz and 5 GHz yielded weak detections. At 5 GHz, the 
pulse was unexpectedly broad, consistent with the results 
of the frequency-subband fit to the pulse shape. Thus, 
in addition to the importance of scatter-broadening at 
the lower frequencies, the pulse appears to be intrinsi- 
cally broad, roughly 30% of the pulse period, full width 
at half maximum. The well-established trend of pulse 
widths increasing as P -0 ' 5 underestimates this result by 
a factor of 2-3 (e.g., for the outermost conal component 
of Mitra k Deshpande 1999). 

Polarization calibration and analysis using the 
PSRCHIVE software package (Hotan et al. 2004) reveals 
that the radio pulse is substantially linearly polarized, 
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despite the depolarizing effects of multipath propagation 
in the scattering tail, incomplete removal of Faraday ro- 
tation (due to the low signal-tonoise ratio), and other ef- 
fects. The true polarized fraction in the pulse is likely to 
be significantly larger than we have observed. The varia- 
tion of the polarization vector’s position angle with rota- 
tional phase, shown in the top-most panel of Figure 2, is 
found to be fairly flat, inconsistent with the canonical 
rotating-vector model (Radhakrishnan & Cooke 1969), 
but similar to behavior seen in the Crab pulsar (espe- 
cially its high-frequency components; Moffett &c Hankins 
1999) and several millisecond-period “recycled” pulsars 
(Ord ei al. 2004). Perhaps not coincidentally, many of 
these same pulsars exhibit nonthermal magnetospheric 
pulsations in X-rays. Finally, we note that the pulsar’s 
rotation measure is RM ~ +270 rad m -2 . 

4 RXTE OBSERVATIONS AND RESULTS 

The field containing CXOU J202221 was observed by 
RXTE for a total of 99 ks over the 8-day span 2010 
January 27-February 4 UT (observation P95316). The 
data were collected with the Proportional Counter Ar- 
ray (PCA; Jahoda et al. 1996) in the GoodXenon mode 
with an average of 1.8 of the 5 proportional counter units 
(PCUs) active. In this mode, photons are time- tagged to 
0.9 /is and have an absolute uncertainty better than 100 
jj , s (Rots et al. 1998). The effective area of five combined 
detectors is approximately 6500 cm 2 at 10 keV with a 
roughly circular field-of-view of ~ 1° FWHM. Spectral 
information is available in the 2-60 keV energy band with 
a resolution of ~ 16% at 6 keV. 

Production PCA data for this observation were ob- 
tained from NASA’s HE AS ARC archive and time-filtered 
using the standard criteria. This rejected intervals of 
South Atlantic Anomaly passages, Earth occultations, 
and other periods of high particle activity, resulting in 
a total of 94.2 ks of useful data. The photon arrival 
times were transformed to the solar system barycenter 
in Barycentric Dynamical Time (TDB) using the JPL 
DE200 ephemeris and the Chandra-derived coordinates 
of CXOU J202221. 

4.1. Timing Analysis 

The timing analysis was restricted to PCA channels 2- 
50, corresponding approximately to the 2-20 keV energy 
range, and from the top Xenon layer of each PCU, to op- 
timize the signal-to-noise for a typical pulsar spectrum. 
Because of the long observation span, photon detection 
times v/ere accumulated in 4-ms bins and a 2 28 -point ac- 
celerated Fast Fourier Transform (FFT) was performed 
to search for a periodic signal, testing a range of plau- 
sible frequency derivatives. This immediately revealed 
a signal at 24 ms of sufficient strength to analyze short 
intervals of data and build up a fullv phase- coherent tim- 
ing solution using the time of arrival (TO A) method, as 
follows. 

The RXTE observations of CXOU J202221 distributed 
over the 8 days were clustered into 10 groups of two, 
three, or four adjacent 96-min orbits, with large gaps in 
between. For each of these groups, we were able to ex- 
tract the period and phase of the signal with sufficiently 
small uncertainties to maintain cycle counts between fit- 
ted epochs. Pulse profiles were generated using the op- 
timum period derived from the Z\ test (Buccheri et al. 
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Figure 3. RXTE timing results for PSR J2022-J-3842 in the 2- 
20 keV X-ray band. Top: Phase residuals after fitting the phase- 
connected quadratic ephemeris given in Table 1 to the set of PCA 
observations as described in §4.1. On average, individual TOA 
measurements were obtained in 17 ks-PCU of exposure. Bot- 
tom: The pulse profile of the RXTE data folded on the best-fit 
ephemeris. Two cycles are shown for clarity. 


1983) with n = 10 harmonics to allow for the sharp pro- 
file. The resulting profiles were cross correlated, shifted, 
and summed to generate a master pulse profile template. 
Individual profiles were then cross-correlated with the 
template to determine the arrival time (time of phase 
zero) and its uncertainty at each epoch. 

The derived TOAs were iteratively fitted to a quadratic 
ephemeris using the TEMPO software. We started by fit- 
ting TOAs from the first three most closely spaced TOAs 
to a linear solution, and then iteratively added the next 
TOA. At each step we found that the new TOA would 
match to < 0.02 cycles the predicted phase derived from 
the previous set. The resulting ephemeris referenced to 
the mid-point of the observation is presented in Table 1 
and the phase residuals are shown in Figure 3. The resid- 
uals are all less than 0.02 cycles and appear to be random 
within their statistical uncertainties, with the exception 
of a single 3a data point (A<£ = 0.025 ± 0.01) which an 
investigation finds no reason to exclude. 

Figure 3 displays the pulse profile using all of the 2- 
20 keV data folded on the final ephemeris. It has a single 
symmetric peak that is triangular in shape and narrow, 
with a FWHM of 0.06 of a full cycle. The measured 
pulsed emission corresponds to 0.91% of the total PCU 
countrate; however, as will be shown below, the intrinsic 
signal is nearly 100% pulsed compared to the flux derived 
from the spectrum of CXOU J202221. We see no energy 
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Tkble 1 

Properties of PSR J 2022 +3842 


Parameter Value 


X-ray Measurements 

R.A. (J2000) a 20 h 22 m 21®689(6) 

Decl. (J2000) a +38°42T4 , /82(7) 

Epoch (MJD TDB) b 55227.00000027 

Period, P (ms) b 24.2877561082(84) 

Period derivative, P h 4.3064(93) X 10” 14 

Range of timing solution (MJD) b 55223-55231 

Radio Measurements 

Long-term average P 4.3192(27) X 10“ 14 

Range of P fit (MJD) 54957-55469 

Dispersion measure, DM (pc-cm -3 ) 429.1 ± 0.5 

Scatter-braodening timescale, (ms) 55(7) 

Flux density, S 2 GHz (/*Jy) 60 

Derived Parameters 

Characteristic age, r c (kyr) 8.9 

Spin-down luminosity, E (ergs -1 ) 1.19 X 10 38 

Surface dipole magnetic field, B 3 (G) 1.03 x 10 12 


Note. — la uncertainties given in parentheses. 
a Chandra ACIS-S3 position registered using 2MASS objects (see 
text). 

b RXTE phase-connected ephemeris. 

dependence of the pulse profile when subdividing the 2- 
20 keV RXTE energy range. 

4.2. RXTE Spectral Analysis 

The pulsed-flux spectrum of PSR J 2022+3842 can 
be isolated through pihase-resolved spectroscopy. We 
used the FTOOLS f asebin software to construct phase- 
dependent spectra based on the ephemeris of Table 1. 
We combined spectra from all observation intervals us- 
ing only data recorded in the top Xenon layer of PCU#2; 
this PCU was (uniquely) used for all observation inter- 
vals over the 8 days. Similarly, we averaged the stan- 
dard PCA responses generated for each epoch. In fitting 
the pulsed flux, the unpulsed emission provides a near 
perfect background estimate. The spectra were accu- 
mulated in two phase intervals representing the off-peak 
(0.1 < <fi 0 ff < 0.88) and on-peak (0.88 < 4> on < 1.1) emis- 
sion and fitted using XSPEC. The fits were performed in 
the 2-16 keV range, above which the background domi- 
nates. 

Based on the Chandra spectrum, we fitted an ab- 
sorbed power-law model with the interstellar absorption 
held fixed at Nh = 1.6 x 10 22 cm -2 ; leaving Nh un- 
constrained results in a larger uncertainty in T. The 
resulting best- fit photon index is T = 1.1 ± 0.2 with 
X 2 = 27.9 for 30 DoF. The absorbed 2-10 keV flux 
for the pulsed emission is 5.4 x 10 -13 erg cm" 2 s _1 , 
which represents all of the point-source flux measured 
from CXOU J202221, indicating that its intrinsic pulsed 
fraction is essentially 100%. In Figure 4, we show the 
result of a joint fit to the RXTE pulsed emission from 
PSR J2022+3842 and CXOU J202221 in the 0.5-16 keV 
band, leaving only their normalizations free. The best 
fit parameters are Nh = (1.7 db 0.3) x 10 22 cm -2 and 
F = 1.0 d= 0.2 with x 2 — 99,3 for 100 DoF. There 
is no indication of spectral curvature in the fitted en- 
ergy range. The independently measured absorbed 2- 
10 keV flux is the same for both CXOU J202221 and the 



Figure 4. The Chandra ACIS spectrum of CXOU J202221 
(black) and the RXTE PCA spectrum of the pulsed emission from 
PSR J2022+3842 (red), fitted jointly to an absorbed power-law 
model, with independent normalizations. The solid line shows the 
best-fit model (see text §4.2). The pulsed RXTE flux is obtained 
by subtracting the off-peak spectrum from the on-peak spectrum. 
Residuals from the best fit are shown in units of standard deviation. 

pulsed emission from PSR J2022+3842, namely Fx = 
5.3 x 10 -13 erg cm -2 s -1 . We adopt the results of the 
joint spectral fit as our final reported values. 

5. DISCUSSION 

There can be little doubt that CXOU J202221 and 
PSR J 2022+3842 are one and the same: the non- thermal 
X-ray source lies at the center of G76.9+1.0, lacks an 
optical counterpart, and anchors nebular emission con- 
sistent with a PWN. Most importantly, the pulsed flux 
detected using RXTE accounts for all of the Chandra 
emission and no more. These properties are consistent 
with pulsars producing broadband magnetospheric emis- 
sion powered by rotational energy losses. Chandra as- 
trometry thus locates the pulsar to sub-arcsec precision. 

PSR J2022+3842 is outstanding in several respects. It 
is the second-most energetic Galactic pulsar known, af- 
ter the Crab pulsar, and the fourth overall, taking into 
account two LMC pulsars, in N157B (PSR J0537-6910, 
with E = 4.9 x 10 38 erg s _1 and P = 16 ms) and in SNR 
0540-69 (PSR J0540— 6919, with E = 1.5 x 10 38 erg s" 1 
and P = 50 ms). However, it is among the least efficient 
at converting its spin-down luminosity into X-rays 1 , with 
r} = L 2 — lokev/^ = 5.9 x 10 -5 D 2 0 . The flat pulsed spec- 
trum also sets it apart from the typical young, energetic 
pulsar, deviating significantly from the observed trend 

T oc 1 /Ve (Gotthelf 2003; 2-10 keV), which predicts an 
index of T = 1.8. Flat power laws have primarily been 
found among pulsars with H.E.S.S. counterparts (e.g., 
PSR J1838— 0655; Gotthelf & Halpern 2008). 

PSR J2022+3842 is also the second-most rapidly rotat- 
ing young pulsar after the LMC pulsar J0537-6910, and 
the shortest-period radio-bright pulsar known. Its ra- 

1 For ease of comparison with the compilation of 
Kargaltsev Sc Pavlov (2008), the flux, luminosity, and efficiency 
in the 0.5-8 keV band are Fpsr — 4.0 x 10 -13 ergs cm -2 o -1 , 
L x = 6.5 x 10 33 Df 0 ergs s -1 , and Lx/E = 5.5 x 10" 5 Df 0 . 
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Figure 5. The true age r of PSR J2022+3842 as a function of 
its spin period at birth, P 0 , according to Equation (1); n is the 
braking index. The curves correspond, from top to bottom, to 
n — 1-5, 2.0, 2.5, 3.0 (dashed), 3.5, and 4.0, converging at the 
present-day period, P = 24.28 ms. The associated upper-limits on 
the pulear age (for Po much smaller than P) are ~ 36, 18, 12, 9, 7, 
and 6 kyr, respectively. 


dio properties are thus potentially interesting probes of 
the elusive radio emission mechanism in pulsars, but de- 
tailed studies will be hampered by this object’s faintness 
and the interstellar propagation effects imposed by its 
great distance. The short rotation period also engenders 
uncertainty in the pulsar’s (and thus the supernova rem- 
nant’s) age: the characteristic age r c = P/2P = 8.9 kyr 
approximates the true age r only when Pq P, where 
Pq is the spin period at birth, 


P 

(n - 1)P 



( 1 ) 


where n is the so-called “braking index.” For a birth pe- 
riod of 16 ms, the shortest known among young pulsars, 
the true age decreases to « 5 kyr, assuming spin-down 
due tc magnetic dipole braking (rc = 3). A counter- 
ing effect, however, is spindown that deviates from the 
standard dipole assumption: the measured braking in- 
dices for several young pulsars are found to lie within 
the range 2.0 < n < 3.0 (see, e.g., Melatos 1997). For 
example, as shown in Figure 5, Pq = 10 ms for rc = 2 
would imply an age of 11 kyr. For any assumed value 
of n, an upper limit to the age of the pulsar, and thus 
the remnant, is provided in the limit P 0 — > 0. A braking 
index well below n = 1.5 would be highly unusual, so 
that ~ 40 kyr is a reasonable upper limit on the age of 
G76.9+1.0. 

PSR J2022+3842 lacks the bright X-ray PWN ex- 
pected from an energetic pulsar: it is only the second ex- 
ample (of some 20) of a pulsar with E £ 4 x 10 36 erg s -1 
unaccompanied by a PWN of comparable brightness, 
Ppwn/Ppsr ~ 1> defying the trend presented in Gotthelf 
(2004) As such, it is similar to PSR J1617— 5055 (Torii 
et al. 1998), the 69 ms pulsar with E = 1.6 x 10 37 erg s~ x , 
B s = 3.1 x 10 12 G, and r c = 8.1 kyr. The un- 
derluminous PWNe in these cases remain unexplained. 
In other respects, however, the X-ray PWN around 
PSR J2022+3842 may not be so unusual. The semi- 
major axis of the elliptical elongation around the unre- 


solved source is approximately 6" in length, representing 
a physical dimension of 9 x 10 17 Dio cm. For compari- 
son, the Vela PWN’s X-ray “outer arc” lies at a distance 
1 x 10 17 cm from the pulsar in the geometric model of 
Helfand et al. (2001). 

If the extended emission around PSR J2022-1-3842 does 
arise from axisymmetric features, the Chandra image 
constrains — in analogy with the Crab, Vela, and others — 
the orientation of the pulsar’s spin axis: in inclination, 
to ~ 50°, and on the sky, to the symmetry axis at po- 
sition angle -64° North through East, orthogonal to 
the long direction of the diffuse X-ray emission and the 
ridge linking the radio lobes. The apparent alignment 
of spin axes and proper motion vectors in some of these 
same pulsars offers a testable prediction for the motion 
of PSR J2022+3842. The spin inclination relative to the 
line of sight, meanwhile, may help elucidate the role of 
viewing geometry (e.g., relative to a narrow magneto- 
spheric emission beam) in i) producing simultaneously 
a narrow X-ray pulse but an unexpectedly broad radio 
pulse, a reversal of their typical relationship, and ii) de- 
termining whether the observed neutron star spectrum 
is primarily thermal or nonthermal: the flat spectrum of 
PSR J2022+3842 departs markedly from the dominant 
thermal emissions of the central stars in the Vela and 
DA 495 SNRs. 

We have already alluded to some of the notable simi- 
larities between G76. 9+1.0 and DA 495, both of which 
can now unambiguously be characterized as radio PWNe. 
Both are very bright radio sources without a well-defined 
outer boundary. Both manifest a bipolar structure and 
have a steep radio spectrum, unusual for a PWN. DA 
495 has a break in its radio spectrum at roughly 1 GHz, 
and G76. 9+1.0 may have one below that as well. Both 
show very faint X-ray emission compared to their ra- 
dio emission, and a small X-ray nebula compared to the 
radio nebula. The conclusions arrived at by ASL+08 
for DA 495 may thus also apply to G76.9+1.0. For 
example, based on the radio/X-ray PWN size ratio of 
DA 495, > 25, and other lines of argument (see also 
Kothes et al. 2008), ASL+08 suggested that the pul- 
sar wind energizing DA 495 has a high magnetization 
factor — i.e., it carries electromagnetic flux and particle 
flux at roughly comparable levels — in contrast with the 
Crab, which has a strongly particle-dominated wind, but 
similar to independent assessments of the wind proper- 
ties of the Vela pulsar. Moreover, ASL+08 speculated 
that, because particle-dominated winds are necessary for 
efficient conversion of wind luminosity to synchrotron lu- 
minosity, PWNe in which Poynting flux is an important 
wind component may be those with dim X-ray PWNe. 

The PWN size ratio for G76. 9+1.0 is ~ 20, comparable 
to DA 495 and nearly two orders of magnitude greater 
than that for the Crab Nebula. (At a distance of 10 
kpc, G76.9+1.0 may be the largest PWN known in our 
Galaxy, with a physical size of 29 x 35 pc, bigger than 
or comparable to MSH 15—57, which is at least 25 pc in 
diameter, Gaensler et al. 2000, and so far believed to be 
the largest Galactic PWN.) If our prior conclusions hold, 
the wind from PSR J2022+3842 should be highly mag- 
netized well beyond its termination shock, contributing 
to the low X-ray conversion efficiency of its PWN. Why 
this might be so is an open question. 

The overarching conclusion of ASL+08 and 
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Kothes et ai. (2008) was that DA 495 is likely to 
be a PWN of advanced age (~ 20 kyr), but having 
evolved without significant interaction with the am- 
bient medium. All of the characteristics we find for 
G76. 94-1.0 similarly indicate a rather old object, yet 
PSR J2022+3842 has all the characteristics of a young 
object: a spin-down age of 20 kyr would require both 
n < 2 and very high spin at birth, P 0 < 5 ms. We do not 
have a ready explanation for this apparent discrepancy, 
which highlights the importance of having uncovered 
the spin and energetic properties of the central pulsar. 

G76. 9+1.0 and its pulsar, J2022+3842, are clearly 
unusual and require further investigation. Additional 
multi- wavelength observations (e.g., an X-ray search for 
the SNR shell, to aid in constraining the system’s age) 
may be fruitful in providing important components for 
their understanding. Based on its spin-down luminos- 
ity, and given the sub-arcsec localization and available 
timing information, PSR J2022+3842 is a good candi- 
date for a search for gamma-ray pulsations using Fermi 
data. The absence of a bright PWN, however, makes it 
an unlikely TeV target, even though the spectrum of the 
pulsed X-ray emission suggests otherwise. 
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